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In 1976 Mandelstam and 't HooftS proposed that the force between a quark- 



ih i ABSTRACT 

0\ ■ 

^ ■ We describe long distance QCD by a dual theory in which the fundamental vari- 

ables are dual potentials coupled to monopole fields and use this dual theory to deter- 
mine the effective Lagrangian for constituent quarks. We find the color field distribution 
surrounding a quark anti-quark pair to first order in their velocities. Using these distri- 

q\ ■ butions we eliminate the dual potentials and obtain an effective interaction Lagrangian 

Lj(xi ,x 2 ',Vi, v 2 ) depending only upon the quark and anti-quark coordinates and ve- 
. locities, valid to second order in their velocities. We propose Lj as the Lagrangian 

describing the long distance interaction of constituent quarks. 

1. Introduction 

m ■ 

m : 
o . 

antiquark pair is analogous to the force between a monopole anti-monopole pair in a 
superconductor. A dual Meissner effect prevents the electric color flux from spreading 
out as the distance between the quark anti-quark pair increases. As a result a linear po- 
tential develops which confines the quarks in hadrons. A concrete model for this "dual 
superconductor" mechanism is realized by a dual theoryi defined by a Lagrangian den- 
• sity C Bfi (C^, Bi) in which the fundamental variables are dual potentials C M coupled to 

>• ! monopole fields denoted by B{. The local coupling of monopole fields to dual potentials 

gives rise to a C M mass and consequently to a dual Meissner effect. The Lagrangian 
density £ cff (C M , BA of the dual theory describes the same long distance physics as the 
original QCD Lagrangian density C QCD defined in terms of Yang Mills potentials 
which are strongly coupled at long distances. 

In this talk we construct C eS , treating quarks as classical point particles having 
coordinates Xi(t) and velocities Vi(t). We solve the resulting classical field equations for 
C M in the presence of a quark anti-quark pair to first order in their velocities. Using 
these solutions we eliminate the dual potentials from £ cff and obtain an effective 
interaction Lagrangian Lj(xi, x 2 ; Vi,v 2 ) = J dxC cS depending upon the quark and anti- 
quark positions and velocities, valid to second order in these velocities. It contains 
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two parameters a s and the string tension a. We use Lj and the canonical formalism 
to construct the quark anti-quark Hamiltonian and obtain a specific quark model for 
calculating the meson spectrum.! 

We begin by showing that dual potentials give an alternate formulation of classical 
electrodynamics. This gives an example where the dual theory and the original Maxwell 
theory describe the same physics at all distances in a situation where the original 
Maxwell potentials and the dual potentials can have different long distance behavior. 

2. Dual Potentials In Electrodynamics 

The inhomogeneous Maxwell equations for a pair of oppositely charged particles 
moving in a relativistic dielectric medium are: 

- -> - 3D 

V-D = p, VxH = J+ — , (1) 

where 

p(x ,t) = e[6(x - x\(t)) - 6(x - x 2 (t))} , (2) 
J(x,t) = e[vi(t)5(x - xi(t)) - v 2 (t)S(x - x 2 (t))] . (3) 
Introduce a line of polarization P s (x ,t) connecting the particles (a Dirac string): 

P s {x,t) = e dy5{x-y{t)), (4) 

Jx 2 {t) 

where the integral dy is along any path L(t) connecting xi(t) and ^(t). (See Fig. 1.) 



Fig. 1. Dirac string connecting oppositely charged particles 
As the particles move, so does the Dirac string. There results a magnetization M s (x ,t), 

M s (x,t) = J edy x -^5{x - y{t)) . (5) 



It is readily checked that 



- - - BP 

p = -V-P s , and J = VxI s + -^. (6) 

at 

Inserting eqs. (6) into (1), we obtain the solution for D and H of the form: 

-, -> -> 8C 

J D = _VxC-P s , H = -VC - — + M S . (7) 

ot 

Gauss' Law and Ampere's Law (1) have thus become kinematic equations whose solu- 
tion, eqs. 7, define dual potentials = (C° , C). 

— * — * — » — * — * 

The constitutive equations B = fiH and E = -D = fiD determine E and B in 
terms of the dual potentials C M and the magnetic permeability \x = K The homogeneous 
Maxwell equations 

dB 

V-B = 0, and Vx E = -— , (8) 

are then dynamical equations determining C M . 

As an example, consider a homogeneous medium with particles at rest. Then the 
equations V x E = 0, E = fiD, and eq. (7) give the following equation for C (denoted 
C dirac = Cd): 

V x (-V xC„)=VxP s . (9) 

Comparing eqs. (9) and (4) with the equation determining the vector potential of a 
magnetic dipole we obtain: 

Jx 2 (t) 4vr \x-yf 
Eqs. (4), (7) and (9) then yield the electric field 

D = -V X C D - P s = ^COULOMB = D c , (11) 

where 



= £ r^V-r^T ■ (12) 

47T | | x — X\ \ \ X — %2 J 



The first term —V x Cd in eq. (11) produces a field analogous to the magnetic field of 
a line of magnetization flowing through the string. The second term — P s cancels the 
field through the string leaving the desired pure Coulomb field. (See Fig. 2.) 

For slowly moving charges, the equations V • B — , B — fiH , V • Cd = and 
eq. (7) give the following equation determining the scalar potential (denoted Cod): 



- V 2 C QD = -V • M s . 



(13) 



Fig. 2. Diagram representing string cancellation mechanism of eq. (11) 



Coo= [:?xfi^. (14) 



Comparing eq. (13) with the equation determining the scalar potential of an electric 
dipole, and using eq. (5), we obtain: 

r*i(*) edy (x - y) 

lx 2 (t) 47r | x — y 

Eqs. (5), (7), (10) and (14) then yield the magnetic field 

-> -> 9Cd -> -> -> 

H = — VC d tt: — I" M s = H BIOT SAVA rt = H BS , (15) 

at 



where 



V\ x (x — Xi) V2 x (a 7 — x 2 ) 



|x-fi| 3 \x-x 2 \ 



(16) 



is the usual Biot Savart magnetic field produced by slowly moving charges. 

To obtain the covariant form of the equations of the dual Maxwell theory we 
define the dual polarization tensor G s a/3 : 

G'lk = M s,k , G\j = -eijkP s ,k ■ (17) 
The field equations (8) can then be written in the form: 

d a ii{d a C p - d p C a ) = -dPfiGlp . (18) 
These equations follow from a dual Lagrangian density Cc'- 

Cc = -^{d a C p -d p C a + G^f. (19) 
Cc describes the same physics as the Maxwell Lagrangian La '■ 

C A = -\{d a A p - d p A a f + J a A a , (20) 

determining the usual Maxwell equations, 

& x e{d a Ap-dpA a ) = J p . (21) 

If e — > at long distances then A^ is strongly coupled (anti-shielding), but since 
fj, — \ — > oo, C M is weakly coupled at long distances. 



3. Non-Abelian Theory 

What is the Lagrangian density C cS describing long-distance QCD in terms of 
dual potentials? Mandelstami has shown that we also have the freedom to use dual 
potentials in non-Abelian gauge theory, 't Hoofti showed that if the Wilson loop 
satisfies an area law then the dual Wilson loop (defined in terms of dual potentials) 
satisfies a perimeter law. This is the non-Abelian analogue of the fact that in a medium 
with e < 1, dual potentials are screened at large distances. Although in non-Abelian 
gauge theory there is no direct relation between ordinary Yang Mills potentials A a 
and dual potentials C°, we can determine the form of L eS by imposing the following 
requirements: 

1. C cS must be invariant under non-Abelian gauge transformations of the Cy 

= Q^C^Q + Ul- l d^VL , (22) 
where = Z)o=i C^X a , Q is an SU (3) matrix and g = 2£ where e is the Yang 

2 

Mills coupling constant, i.e. a s = |-. 

2. C cS generates a mass for the field via a Higgs mechanism coupling C M to 3 
scalar octets Bi(i = 1,2,3) carrying monopole charge. 

3. We write a minimal C eS satisfying (1) and (2). 

4. We assume higher dimension operators not included in £ cff are not quantitatively 
important at large distances. 

First consider the coupling of a quark anti-quark pair qq to C M . We choose a 
gauge where 

p(x) = eY[5(x — xi) — 8(x — X2)) , (23) 

where Y is the hypercharge matrix having diagonal elements 3 , | , — | and vanishing 
off diagonal elements (2trF 2 = |). Then the covariant polarization tensor G s Q/3 for 
the Dirac string connecting the quark anti-quark pair is given by eqs. (4), (5), and 
(17), multiplied by Y. The dual potentials are then also proportional to Y and so 
[C M , C u ] = 0. One unit of gauge invariant Z 3 flux flows from q to q along the Dirac 
string. To see this, evaluate 

= p e t9 $ &<a ^ (24) 

where we choose the path in (24) to be a small circle surrounding the Dirac string. 
Near the string the effect of the coupling to the monopole fields Bi can be neglected 
and C ~ YCd- Then eq. (24) becomes, 



(25) 



Fig. 3. Quark anti-quark pair moving along trajectories x\{t) and x 2 (t), respectively. 

— * 

Any continuous deformation of C in SU (3) leaves $ unchanged. 
We choose a gauge for the monopole fields Bi so that, 

B 1 = A 7 B 1 {x), B 2 = -\ 5 B 2 (x), B 3 = \ 2 B 3 (x). (26) 

As x — > oo, Bi(x) — > So where So is the vacuum value of the Bi determined by the 
position of the minimum of the Higgs Potential W(B): 

5W 

— = 0, a tB l = B . (27) 

ODi 

There is no SU(3) gauge transformation which leaves invariant the vacuum values of all 
three B^. Dual SU(3) symmetry is then completely broken and all the dual potentials 
become massive. 

The explicit form of the dual Lagrangian C cfi coupling a quark anti-quark pair to 
is then given by, 

C cH = 2 tr {-1 (0„C„ - d u C, + G%f + liV.B,) 2 } - W , (28) 

where 

V^B^d^-igiC^Bi]. (29) 

Using eqs. (26) and the fact that = YC^(x) we can eliminate the color matrices 
from eq. (28) to obtain 

C ee = - d u C, + G%f + 2[(d fl B 1 ) 2 + (d,B 2 ) 2 + (d,B 3 ) 2 + g 2 C\B\ + B\)\ - W . 

(30) 

4. Effective Lagrangian for Constituent Quarks 

Consider a quark (anti-quark) at Xi(t)(x 2 (t)) moving with velocities Vi(t)(v 2 (t)) 
(see Fig. 3). To find the effective interaction Lagrangian Li(x\ , x 2 ; vi , v 2 ) we solve the 
classical field equations for C M and Bi generated by £ eff , eq. (30), to first order in V\ 
and v 2 , and obtain the classical solutionsJl 

C M = C^x,R(t) , v^t) , v 2 (t)) = , B % = Bi(x,R(t) , v^t) , v 2 (t)) = Bf- , 

(31) 



Fig. 4. Static flux tube configuration 



where R = x\ — x 2 . 

For v\ — v 2 — 0, the scalar potential C*Q lass = 0, and C cl ^ s and Bf^ reduce to 
static flux tube solutions C st3lUc (x , R), £?f tatlc (x , R), describing a configuration where 
the color flux! connecting the quark anti-quark pair is confined to a tube of radius Rft 
(see Fig. 4). The monopole fields vanish at the center of the flux tube and approach 
their vacuum value B at large distances. We insert this static solution into £ cff and 
integrate over all space to obtain the central potential Vq(R): 

V (R) = - f dxC cS (C = C static , C = , Bi = 5f atic ) . (32) 

At large R, Vq(R) — > crR, where the calculated string tension a ~ (5-Bo) 2 - The experi- 
mental string tension a ~ .2GeV 2 then determines Bq. 

To first order in v\(t) and ^(i), the static distributions, (7 static (f , R) and Bf at[c (x , R) 
follow the quark motion adiabatically. The time dependence of (7 class and Bf^ thus 
arises only from the explicit time dependence of R. The resulting time varying color 
electric fields along with the moving Dirac string then generate a scalar potential Co lass 
and color magnetic fields. These produce a velocity dependent term — V 2 in the full 
quark anti-quark effective Lagrangian Lj given by: 

L T {x ! , x 2 ;vi , v 2 ) = J dxC cS {C^ = C^ ss , B { = Bf*«°) , (33) 

valid to second order in vi and v 2 since the Lagrangian is stationary about solutions 
to the static field equations. 

Separating Lj into its color electric and color magnetic components, we have 

-L I (R,v 1 ,v 2 ) = V (R) + V 2 (R,v 1 ,v 2 ), (34) 

where the static potential Vo(-R) is determined by eq. (32). The color magnetic com- 
ponent V 2 , obtained from eq. (33), has the form! 

v,(jU,%) = [Sx f R : vW V4R)+ lS - (S ^ )] \ ,(R) (35) 



The spin dependent potential V^ P i n includes both spin-spin and spin orbit contri- 
butions, i The spin orbit contribution^ has the qualitative features of Monte Carlo data 
obtained from lattice gauge theory.0 The potentials V-(R) , V\\(R) , V + (R) , and Vl(R) 
are moments of the field distributions determined by the integrand in eq. (33). The 
potentials V-(R) and V[i(i2) are related to Vq(R) by Lorentz invariance:! 

V.(K) = ~V (R) , V \\(R) = ~\vo(R) + ■ (36) 

The function Vl(R) determines the energy associated with longitudinal oscillations 
R • The function V + (R) determines the moment of inertia I(R) of the rotating flux 
tube according to the equation: 

I(R) = -\r 2 V + {R) . (37) 

At large R, V+(R) — > —AR, where A ~ .21a. This result for V + (R) along with the 
large R result V (R) — ► aR gives a leading Regge trajectory! a(M 2 ) which is linear for 
large M 2 with a slope a « « .75GeV~ 2 . 

We list below analytic parameterizations of the integrals over the field distribu- 
tions determining V (R) , V+(R) and Vl(R): 

Vo(R) = -^e- 51 ^ + aR - Q.Qh{a s a) 1 ' 2 , (38) 

o JrC 

V+(R) = -^e" 1 ' 13 ^ - .21aR + 1.12^^, (39) 

6K 

Vl{R) = -^e- 68 V^ R + .09^ . (40) 
6K 

These potentials determine generalized Wilson loops with additional electric and mag- 
netic field insertions and hence can be compared with Monte Carlo calculations of these 
quantities once they are carried out. 

For small R the effective Lagrangian Lj (eqs. (34)- (40)) approaches Ljj, the Dar- 
win Lagrangian (multiplied by the color factor |) describing the interaction of electrons 
and positrons after the elimination of the electromagnetic field.0 Thus in the small R 
limit Li gives the one gluon exchange potential to order {^-) 2 just as Darwin's classical 
calculation gives the one photon exchange potential to the same order. Lj can therefore 
be extrapolated smoothly from the large R confinement region to the short distance 
perturbative domain. It cannot be used at shorter distances where radiative corrections 
giving rise to asymptotic freedom and a running coupling constant must be accounted 
for. 

It should be emphasized that the potentials (38)- (40) are obtained from the inte- 
grated field energy associated with the color field distributions surrounding the moving 



quark anti-quark pair. As R increases, the field lines are compressed and the result- 
ing Lj decreases less rapidly then Lj> At large R the color fields evolve into a flux 
tube distribution and give rise to the terms in Lj linear in R. Lj is not calculated 
as a superposition of a long distance confinement contribution and a short distance 
perturbative term as eqs. (38)-(40) might suggest. These formulae are just analytic 
parameterizations of the numerical integrals described above and the R dependence of 
these potentials reflect the evolution of the color field distributions with R. 

Finally we note that the whole concept of an effective Lagrangian depending 
only upon particle positions and velocities has a meaning only when radiation can 
be neglected. To estimate when radiation is important, consider a quark anti-quark 
pair separated by a distance R and moving on a circular orbit with frequency uj. 
Then radiation occurs if tkv > Mc 2 , where M is the mass either of the C M field or of 
the monopole field. (These masses are of order 500 MeV). Since u ~ j|, this means 
he/ R > Mc 2 , or R < Ti/Mc. Thus, for R > -p- ~ Rft, radiation should be suppressed 
and the quark anti-quark effective Lagrangian can be used. 

5. Energy Levels of bb and cc Systems 

Starting with Lj augmented by non-relativistic quark kinetic energy terms, we 
have constructed the Hamiltonian H for a heavy quark anti-quark system by the canon- 
ical procedure. We determined the parameters a s , a , m c and m h by finding a best fit! 
to the 17 known levels of cc and bb systems. Our best fit parameters are a s = .37, 
cr = .20GeV 2 , m c = lSAGeV, nib = 4.77GeV. These values of a s and a lead to a flux 
tube radius Rft ~ .5fm. 

6. Summary 

(i.) We have proposed that the long distance physics of QCD is determined by a 
dual theory describing the interactions of dual potentials and monopole fields Bi 
which are weakly coupled at long distances. The physical equivalence at long dis- 
tances of £ cff (C M , Bi) and C QCD depending upon strongly interacting allows us to 
use C cS (Cn, Bi) to calculate the long distance properties of QCD (electric magnetic 
duality) B 

(ii.) We have used C eS to construct an effective Lagrangian Lj for constituent quarks, 
which leads to a definite quark model that can both be compared with experiment and 
with lattice gauge theory simulations. 
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